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Lipid liposomes are promising drug delivery systems because they have superior 
curative effects owing to their high adaptability to a living body. Lipid liposomes 
encapsulating proteins were constructed and the structures examined using 
synchrotron radiation small- and wide-angle X-ray scattering (SR-SWAXS). 
The liposomes were prepared by a sequential combination of natural swelling, 
ultrasonic dispersion, freeze-throw, extrusion and spin-filtration. The liposomes 
were composed of acidic glycosphingolipid (ganglioside), cholesterol and 
phospholipids. By using shell-modeling methods, the asymmetric bilayer 
structure of the liposome and the encapsulation efficiency of proteins were 
determined. As well as other analytical techniques, SR-SWAXS and shell- 
modeling methods are shown to be a powerful tool for characterizing in situ 
structures of lipid liposomes as an important candidate of drug delivery systems. 
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1. Introduction 

Liposomes are closed bilayer lipid systems and afford a very 
useful tool in various scientific fields including not only 
biology but also theoretical physics, biophysics, chemistry, 
colloid science and so on. In addition, during the past 30 years, 
liposomes have received a lot of attention as effective drug 
delivery systems (DDS) because they can reduce drug toxicity 
due to biodegradability and biocompatibility, and offer 
promising carriers of anti-cancer, anti-fungal and anti-biotic 
drugs, gene medicines and anesthetics and anti-inflammatory 
drugs, compared with other delivery systems (Kaneda, 2000; 
Rahimpour & Hamishehkar, 2012; Allen & Cullis, 2013). 
Liposome studies and those clinical trials have progressed 
from the use of conventional liposomes to that of 'second- 
generation liposomes' that are developed by modulating the 
lipid composition, size and modified surface of liposome 
(Torchilin, 2005). As a candidate reagent of surface modifi- 
cation of liposomes, glycolipids or sialic acids are considered 
to be promising. 

On the other hand, the formation of a lipid microdomain in 
mammalian plasma membrane, a so-called lipid raft (Simons 
& Ikonen, 1997, 2000; Hakomori, 2001; Simons & Toomre, 
2001), has been attracting intensive interest because lipid rafts 
are assumed to have functions as platforms of membrane- 
associated events such as signal transduction, cell adhesion, 



lipid/protein sorting and so on (Anderson & Jacobson, 2002; 
Gassart et al, 2003; Helms & Zurzolo, 2004). A common 
feature of lipid rafts is their peculiar lipid composition, being 
rich in glycosphingolipids (GSLs), sphingomyelin and 
cholesterol. Gangliosides, major components of GSLs, are 
acidic lipids composed of a ceramide linked to an oligo- 
saccharide chain containing one or more sialic acid residues, 
which are abundant in central nervous systems. Functions of 
lipid rafts are assumed to relate closely to the peculiar features 
of GSL molecules both in ceramide and oligosaccharide 
portions that can form complex hydrogen-bonding networks 
(hydrogen-bond donor and acceptor) (Pascher, 1976). 
Therefore, it can be assumed that gangliosides are promising 
reagents of liposome due to those intrinsic properties. 

By using neutron and synchrotron X-ray scattering techni- 
ques [small-angle X-ray scattering (SAXS), small- angle 
neutron scattering (SANS)], we studied various functional 
properties of gangliosides and those aggregates with other 
lipids under various conditions. The temperature dependence 
of the hydration of the sugar head region and the dissociation 
degree of sialic acids (Hirai et al., 1996a,5,c, 1999; Hirai & 
Takizawa, 1998; Hayakawa & Hirai, 2002), the maximum 
miscibility of cholesterol molecules against gangliosides and 
the cholesterol-dependent micelle-to-vesicle transition (Hirai 
et al., 2005), the vesicle-to-lamellar reversible transition by 
Ca 2+ ions (Hayakawa & Hirai, 2003), and the asymmetric 
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bilayer structures of ganglioside/cholesterol/phospholipid 
liposomes (Hirai et al, 2003; 2006) were clarified. In addition, 
by using the neutron spin-echo technique, we found the 
function of ganglioside molecules controlling the undulation 
motion of micelles (Hirai et al., 2001) and the bending 
modulus of liposomes (Hirai et al., 2005). Recently, we have 
found that the interaction between ganglioside and amyloid 
beta protein (A^l-40) significantly reduces the diffusional 
bending motion of liposomes (Hirai et al., 2013). These 
previous results show that lipid-rafts' components (ganglio- 
side/cholesterol-rich microdomains) can modulate properties 
of membrane structure and dynamics through the interaction 
between lipid and protein (hydrophobic coupling between a 
protein and a surrounding lipid bilayer) (Lundbaek et al., 
2004), suggesting that ganglioside molecules are useful for 
designing a new type of liposome DDS. 

Based on the above results, we have studied the structural 
characteristics and stability of lipid-rafts model liposomes 
entrapping proteins by using synchrotron radiation small- 
angle and wide-angle X-ray scattering (SR-SWAXS). The 
liposomes investigated in the present study would be 
promising components of DDS because they are easily 
adjusted to attaining appropriate curative effects by modu- 
lating lipid composition, size and surface charge of liposomes 
due to their high and specific adaptability to a living body. 
In addition, in spite of numerous studies of liposomes, char- 
acterization of liposomes entrapping proteins, as far as we 
know, has been rarely carried out using scattering methods. 



2. Materials and methods 
2.1. Samples 

The liposomes used were lipid mixtures composed of 
disialoganglioside (G D1 ) from bovine brains, cholesterol and 
phospholipids (egg-PC). Cholesterol was purchased from 
Sigma Chemical Co. (USA). G D1 and egg-PC were purchased 
from Avanti Polar Lipids Inc. (USA). The protein was 
myoglobin from horse skeletal muscle (Sigma). Liposome 
samples were prepared by the sequential combination of 
natural swelling, ultrasonic dispersion, freeze-throw, extrusion 
and spin-filtration, as described below. As a first step, G D1 , 
phospholipid lipids and cholesterol were separately dissolved 
in a chloroform/methanol mixture solvent [1/1 (v/v)]. These 
solutions were mixed in the molar ratio [G D1 ]/[cholesterol]/ 
[phospholipid] 0.1/0.1/1. After mixing, to remove the organic 
solvent, the lipid mixture solutions were dried under a 
nitrogen stream and annealed in vacuo overnight at 318 K. 
The dried mixtures were suspended in myoglobin solutions 
(1%, 2.5%, 5% w/v myoglobin in 10 mM Hepes buffer at pH 
7.4), where the phospholipid concentration was 2% w/v. Those 
suspensions were stored all night to obtain giant multilamellar 
vesicles (GMLVs) by natural swelling. As a second step, the 
suspensions of GMLVs were sonicated for 10 min at ~316 K 
to obtain homogeneous small unilamellar vesicles (SUVs) by 
using a high-power probe-type ultrasonicator (Model UH-50, 
SMT Co.). As a third step, these SUV solutions were subjected 



to a ten-time freezing-thawing process to obtain giant uni- 
lamellar vesicles (GUVs). As a fourth step, large unilamellar 
vesicle (LUV) solutions were prepared by an extrusion 
method using the LiposoFast Basic extruder system (Avestin, 
Canada) with a polycarbonate filter (pore diameter 100 nm). 
The samples were subjected to ~50 passes through the filter. 
As a final step, the obtained LUVs solutions were subjected to 
spin-filtration for removing un-entrapped proteins and for 
concentrating LUVs by using an ultra-filtration device 
(Vivaspin-50K by Sartorius Co.) at 1000 r.p.m. for 10 min x 
5 times. The validity of each step was checked by X-ray scat- 
tering measurement to optimize the above sample preparation 
condition. 



2.2. X-ray scattering measurements 

SR-SWAXS experiments were performed using the BL- 
40B2 spectrometer at the Japan Synchrotron Radiation 
Research Institute (JASRI, Harima, Japan) and by using the 
BL-10C spectrometer at the High Energy Accelerator 
Research Organization (KEK, Tsukuba, Japan). The X-ray 
wavelengths and the sample-to-detector distances were 0.75 A 
for a 51 cm camera length and 1.0 A for a 4089 cm camera 
length at BL-40B2, and 1.49 A for a 190 cm camera length at 
BL-10C. X-ray scattering intensity was recorded using the R- 
AXIS IV IP detector from Rigaku. The exposure time was 10 s 
at BL-40B2 and 180 s at BL-10C. Each sample solution was 
contained in a cell composed of a pair of thin quartz windows 
with 1 mm path length and was served for SWAXS measure- 
ments. 



2.3. Scattering data treatment and shell-modeling analysis 

The background correction of wide-angle scattering data 
was described elsewhere (Hirai et al., 2002, 2004). The distance 
distribution function, p(r), was obtained by Fourier transform 
of the observed scattering intensity, I(q), as 

J oo 

p(r) = ^f rql{q) sin(rq) dq, (1) 

where q = (47r/A)sin(#/2), 0 is the scattering angle and A is the 
X-ray wavelength. We used the following model scattering 
function, I(q,R), of a multi-shelled ellipsoidal particle with a 
size distribution as shown previously (Hirai et al., 2003), 

oo 

Kq)= f I s (q,R)D(R)dR, (2) 

^min 

where R min is the lower limit of the particle radius determined 
by the bilayer thickness of the LUV; D(R) is the number 
distribution function of the particle radius R; I s (q,R) is the 
spherical averaged scattering function of an ellipsoidal particle 
with radius R composed of n shells [ith shell with average 
excess scattering density p i (so-called contrast), radius R t and 
Vi]. I s (q,R) and D(R) are given as follows, 
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I s (q,R)= }h\hVih(qRi)/(aRi) 

0 1 L 

+ t(Pi - Pi-dVihiqRd/iqRi)]] 2 dx 9 (3) 

i = 2 JJ 

where j\ is the spherical Bessel function of the first rank. R t is 
defined by 

R i = ri [l +x ^-l)f\ (4) 

where r t and v f are the semi- axis and its ratio of the zth ellip- 
soidal shell, respectively. For a spherical-shelled particle (v t = 
1, R t ■ = r f ), equation (3) is simplified to 

/ s (^,i?) = 9[p l y l y 1 (^ 1 )/(^i) 

+ E(A " A_i)KA(^)/(^)f. (5) 

i = 2 J 

As a function of D(R), we adopted the Gaussian distribution 
function given by 

D ( R ) = tt-Jtt^pR* - R) 2 / 2 * 2 ], (6) 

(Z7T) CT L -I 

where R and cr are the average radius and the standard 
deviation, respectively. 

3. Results and discussion 

3.1. Removal of un-entrapped proteins from LUV solution by 
spin-filtration and its effect on LUV structure 

In the present sample preparation, not only occlusion of 
proteins within LUVs but also removal of un-entrapped 
proteins from LUV solutions is important. Fig. 1 shows the 
scattering curves of the mixtures of LUV (un-entrapping 
proteins) solution and protein solution before (Fig. la) and 
after (Fig. lb) spin-filtration. The molar ratio of lipids is [G D1 ]/ 
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Figure 1 

Scattering curves of the mixtures of LUV (un-entrapping proteins, [G m ]/ 
[cholesterol]/[egg-PC] = 0.1/0.1/1) solution and protein solution at 298 K 
before and after spin-filtration for removing proteins existing in the 
external solution of LUV. (a) and (b) correspond to before and after spin- 
filtration, respectively. Before spin-filtration in (a), the protein concen- 
trations were 0, 0.5, 1.25 and 2.5% (w/v), the LUV concentration 1.0% 
(hVv). The scattering curves are shifted along the vertical axis to avoid 
overlap. 
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Figure 2 

Scattering curves of LUVs entrapping proteins ([G D1 ]/[cholesterol]/[egg- 
PC] = 0.1/0.1/1) at 298 K before (a) and after (b) spin-filtration. The 
charged (initial) protein concentration was varied from 1% (w/v) to 5% 
(w/v), the LUV concentration 2.0% (w/v). The scattering curves are 
shifted along the vertical axis to avoid overlap. 

[cholesterol]/[egg-PC] = 0.1/0.1/1. In Fig. 1(a), the concentra- 
tion of the LUV was fixed to be 1.0% (w/v), and that of the 
protein was varied from 0.5 to 2.5% (w/v). As shown in 
Fig. 1(a), the presence of amounts of un-entrapped proteins in 
the external solution significantly changes the scattering curve 
in the q ranges of - 0.01-0.15 A -1 and of ~l-2 A" 1 . On the 
other hand, in spite of the addition of proteins with different 
concentrations, all scattering curves after spin-filtration agree 
with the scattering curve of LUV without the addition of 
proteins in the whole q range, indicating that un-entrapped 
(free) proteins were removed from the solution without 
affecting the LUV structure and its dispersity. Figs. 2(a) and 
2(b) show the scattering curves of LUVs ([G D1 ]/[cholesterol]/ 
[egg-PC] = 0.1/0.1/1) entrapped proteins before and after spin- 
filtration, respectively. The charged protein concentration was 
varied from 1% (w/v) to 5% (w/v). After spin-filtration, all 
scattering curves with different protein concentrations show a 
similar profile, suggesting that the proteins un-entrapped in 
LUVs were removed from the solutions. Hereinafter, we call 
LUV un-entrapped proteins and LUVs entrapped proteins 
empty-LUV and filled-LUV, respectively. 

3.2. Structural characteristics of LUVs entrapping proteins 

Fig. 3 shows the scattering curves of the empty-LUV and the 
filled-LUV. The difference between the scattering curves 
appears in the q range of ~ 0.01-0.15 A -1 , which is attribu- 
table to the change of the average excess scattering density 
(so-called contrast) of the core (water pool region) of the 
LUV by the protein occlusion, as shown in the following 
section. The effect of protein occlusion on the LUV structure 
is also seen in the distance distribution function, p(r), obtained 
by using equation (1). Fig. 4 shows the p(r) functions in the 
empty-LUV and the filled-LUV. The position of the maximum 
of the broad peak shifts from ~715 A for the empty-LUV 
to ~590A for the filled-LUV, suggesting that the average 
contrast inside of the LUV becomes higher. In addition, there 
exists another evident difference in the p(r) functions at the 
short-distance region as shown in the insert in Fig. 4. The 
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Figure 3 

Comparison of the scattering curve of LUV un-entrapping proteins 
(empty-LUV) with that of entrapping ones (filled-LUV) at 298 K from 
Fig. 2 {[G D1 ]/[cholesterol]/[egg-PC] = 0.1/0.1/1, charged protein concen- 
tration was 5% (w/v)}. 
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Figure 5 

Experimental and theoretical scattering curves of the empty-LUV ([G D1 ]/ 
[cholesterol]/[egg-PC] = 0.1/0.1/1) at 298 K using equations (2) and (5), 
where the insert is the obtained size distribution function. The reliability 
factor R is 0.051. 



oscillation profile for the filled-LUV in Fig. 4 below ~ 55 A 
tends to be damped compared with that for the empty-LUV. 
The above changes in the scattering curves and in the p(r) 
functions are reasonably explained by the protein occlusion, 
as shown below. 

3.3. Shell-modeling analysis and simulation of LUVs 

Owing to the simple structural symmetry of vesicles, the 
shell-model scattering function given by equations (2) and (5) 
is applicable for describing experimental scattering curves in 
many cases (Hirai et al, 2003, 2013; Hirai, 2007; Onai & Hirai, 
2010). Fig. 5 shows the experimental and theoretical scattering 
curves of the empty-LUV ([G D1 ]/[cholesterol]/[egg-PC] = 0.1/ 
0.1/1), where the insert is the obtained size distribution func- 
tion. The reliability factor R defined by R = E|/ e xperiment(g) — 
Aheory(<3 r )|/^4x P eriment(<2 r ) is 0.051, which indicates the validity 
of the shell-model fitting. The deviation between the theore- 

0.101 1 1 1 1 1 




0 500 1000 1500 2000 
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Figure 4 

Comparison of the distance distribution function, p(r), of the empty-LUV 
with that of the filled-LUV obtained from Fig. 3. The insert expands the 
p(r) functions in the short-distance region. 



tical and experimental S AXS curves above q = ~ 0.2 A -1 is 
attributable to the simplification of the spherical shell struc- 
ture. Fig. 6 shows the contrast profile of the intra-bilayer 
structure of the empty-LUV obtained from the modeling 
analysis. The low positive-contrast region at the outer-leaflet 
of the bilayer corresponds to the protruded portion of the 
oligosaccharide chains from the membrane surface, which is 
the common feature for LUVs containing ganglioside mole- 
cules (Hirai et aL, 2003, 2006, 2013; Hirai, 2007; Onai & Hirai, 
2010). Such an asymmetry of the bilayer structure results from 
the difference in the geometrical packing parameters of lipid 
molecules. Namely, owing to the presence of the oligo- 
saccharide chain containing sialic acid residues in the hydro- 
philic polar head region, gangliosides are preferentially 
localized at the outer-leaflet of the bilayer membrane to 
minimize the surface free energy. The high positive-contrast 
regions reflect the phosphocholine head portions at both sides. 
The positive-contrast region at the inner-leaflet of the bilayer 
corresponds to a high-density hydration shell around the 
phosphocholine heads. 

Based on the above structural parameters of the open-LUV, 
we can simulate the effect of protein occlusion on the scat- 
tering curve of the LUV by using equations (1) and (5). The 
structure parameters of the bilayer were as in Fig. 6. The 
average radius and the standard deviation of the radius were 
400 A and 30%, respectively. The average scattering density of 
water for X-rays is 9.38 x 10 10 cm 2 , and those of proteins 
range from ~11.7 to ~12.0 x 10 10 cm 2 (Stuhrmann & Miller, 
1978). Therefore, the occlusion of proteins increases the 
contrast of the water pool of the LUV from zero to a positive 
value. Fig. 7 shows the shell-model scattering function of LUV 
with a size distribution depending on the change of the water- 
pool contrast accompanied by the protein occlusion. In Fig. 7 
the relative value of the water-pool contrast was varied from 0 
to 2.6, which corresponds to the encapsulation efficiency of the 
protein within the water-pool of LUV from 0 to ~100% (v/v). 
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Figure 6 

Structure of the bilayer of the empty-LUV obtained from the shell- 
modeling analysis in Fig. 5. The obtained value of contrast in each region 
is on a relative scale. 



As shown previously (Onai & Hirai, 2010), the scattering 
curve of an empty-LUV can be characterized by the shoulder 
at ~0.01 A -1 and the broad rounded peak at ~ 0.04-0.2 A" 1 . 
The former and the latter correspond to the spherical shape 
of LUVs with a size distribution (polydispersity) and to an 
internal structure of the bilayer, respectively. With ascending 
encapsulation efficiency, the scattering curve changes 
systematically, especially in the q range below ~0.1 A -1 . The 
shoulder at q = ~ 0.01 A -1 and the dip at q = ~ 0.04 A -1 
become smeared and/or shallower. These changes of the 
scattering curve can qualitatively explain the difference shown 
in Fig. 3. Fig. 8 shows the p(r) functions obtained from Fig. 7. 
With the rise of the water-pool contrast from 0 to 2.6 [the 
encapsulation efficiency of the protein from 0 to ~100% 
(v/v)], the position of the peak maximum significantly shifts 
from ~700 A to ~500 A, and the ripple profile below ~55 A 
becomes smeared out as shown in the insert of Fig. 8. The 
above results of the shell-modeling simulation qualitatively 
describe the tendency of the observed changes in Figs. 3 and 4. 




Figure 7 

Shell-model simulation of LUV scattering function depending on the 
protein occlusion. The encapsulation efficiency of proteins [0 to ~100% 
(v/v)] within the water pool of LUV corresponds to the change of the 
relative contrast of the innermost core of LUV from 0 to 2.6. The change 
of the scattering curve can explain the difference shown in Fig. 3. 
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Figure 8 

Simulated distance distribution function, p(r), of LUV depending on the 
protein occlusion, obtained from Fig. 7. The insert expands the p(r) 
functions in the short-distance region. The changes in the broad peak 
position at ~ 500-700 A and in the ripple profile below ~ 55 A can 
describe the experimental ones in Fig. 4. 



3.4. Estimation of encapsulation efficiency of proteins 

Radii of gyration obtained from small-angle X-ray scat- 
tering data using Guinier plots (Guinier & Fournet, 1955) are 
well known to be useful for characterizing structures of 
particles with random orientation. Fig. 9 shows the water-pool 
contrast dependence of the radius of gyration, R gJ of LUV 
obtained from Fig. 7. In Fig. 9, the R g of LUV decreases 
sensitively against the rise of the contrast (encapsulation 
efficiency), especially at low protein occupancy. In the present 
sample preparation, we varied the charged amount of proteins 
within the LUV from 1% (w/v) to 5% (w/v). Fig. 10 shows the 
experimental R g values of the LUVs obtained from Fig. 2(b) 
for the different charged amount of proteins. The R g value 
decreased from 587 ± 3 A for the empty-LUV to 506 ± 9 A 
[1% (w/v)], to 509 ±7 A [2.5% (w/v)] and to 501 ±6 A [5.0% 
(w/v)] for the filled-LUVs. The experimental R g values of the 
filled-LUVs are similar within experimental error. The insert 
of Fig. 10 shows the simulated relation between the encapsu- 




0.5 1 1.5 2 

change of contrast (arb. units) 

Figure 9 

Radius of gyration, R g , of LUV depending on the protein occlusion, 
which was obtained from the simulated scattering functions in Fig. 8. The 
R g value is plotted against the change of the relative contrast of the 
innermost core of the LUV. 
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Figure 10 

Experimental radii of gyration, R g , of filled-LUVs for the different 
charged amounts of proteins within the LUVs from 1% (w/v) to 5% 
(w/v). The insert shows the relation between the encapsulated protein 
concentration [% (w/v)] and the decrement rate of R g obtained from 
Fig. 8, which was used to estimate the encapsulated protein concentra- 
tion. 

lated protein concentration [% (w/v)] and the decrement rate 
of R g of filled-LUV normalized by R g of the empty-LUV. From 
this relation, we can determine the encapsulated protein 
concentrations of the filled-LUVs to be ~5.5% (w/v). In spite 
of the difference of the initial protein concentration, the 
encapsulation efficiency tends to be higher than the charged 
amount of proteins, suggesting that the present LUV 
preparation induces a slight condensation of proteins within 
the water-pool of the LUV. 



4. Conclusion 

By the combination of different types of liposome preparation 
methods, proteins were effectively encapsulated within the 
water pool of LUV composed of glycosphingolipid (ganglio- 
side G D1 ), cholesterol and phospholipid. Ganglioside mole- 
cules would be promising reagents for designing a new type of 
liposome DDS since gangliosides can modulate the size and 
surface charge of liposomes owing to these intrinsic properties. 
The observed changes in the SWAXS curve of the LUV 
caused by the encapsulation of proteins were reasonably 
explained by the shell-model simulation. The encapsulation 
efficiency of proteins within LUV was also determined, 
suggesting that the concentration of the encapsulated proteins 
tends to be higher than that of the initially charged proteins 
(condensation effect). This would result from the LUV 
preparation method used here. On the other hand, recently, 
different analytical techniques, such as scanning probe 
microscope, fluorescent probe microscopy and so on (Spyr- 
atou et al, 2009), can be applied to characterize liposomes. As 
well as other techniques, the present results show that struc- 
tures of LUVs of lipid mixtures as an important candidate of 
DDSs can be characterized well by using SR-SWAXS and 
shell-modeling methods. 



The X-ray scattering experiments were performed under 
the approval of the program advisory committee of the Japan 
Synchrotron Radiation Research Institute (proposal Nos. 
2009B1730 and 2011A1929) and under the approval of the PF 
program advisory committee (proposal Nos. 2009G202 and 
2010G056). This research project was supported by the 
Grants-in-Aid for Scientific Research of JSPS (the Japan 
Society of the Promotion of Science) (proposal No. 22570152). 
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